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Electric field produced inside a solute by a uniformly polarized liquid is strongly affected by dipolar polar- 
ization of the liquid at the interface. We show, by numerical simulations, that the electric "cavity" field 
inside a hydrated non-polar solute does not follow the predictions of standard Maxwell's electrostatics of 
dielectrics. Instead, the field inside the solute tends, with increasing solute size, to the limit predicted by the 
Lorentz virtual cavity. The standard paradigm fails because of its reliance on the surface charge density at 
the dielectric interface determined by the boundary conditions of the Maxwell dielectric. The interface of a 
polar liquid instead carries a preferential in-plane orientation of the surface dipoles thus producing virtually 
no surface charge. The resulting boundary conditions for electrostatic problems differ from the traditional 
recipes, affecting the microscopic and macroscopic fields based on them. We show that relatively small differ- 
ences in cavity fields propagate into significant differences in the dielectric constant of an ideal mixture. The 
slope of the dielectric increment of the mixture versus the solute concentration depends strongly on which 
polarization scenario at the interface is realized. A much steeper slope found in the case of Lorentz interfacial 
polarization also implies a higher free energy penalty for polarizing such mixtures. 
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I. INTRODUCTION 

When an interface is created in a dielectric, the surface 
dipoles change their preferential orientations relative to 
the dipoles in the bulk. The response of the dielectric 
to a weak external field is then a composite result of the 
response of these surface dipoles and the bulk dipoles. 
The question of whether the dielectric response of a ma- 
terial is sensitive to its surface structure is decided by the 
relative weights of these two contributions. 

The Maxwell electrostatics of dielectrics neglects the 
structure of the interface and replaces the interfacial re- 
gion of a finite microscopic dimension with an infinites- 
imally thin mathematical surface. This mathematical 
surface cuts through the dipoles of the medium creat- 
ing a surface charge with the charge density op (Fig. 
UK). It is given by the projection of the dipolar polariza- 
tion at the interface P(rg) on the outward normal n to 
the surface bounding the dielectric^ o\p(rs) = P„(rs), 
P n (v s ) =n-P(r s ). 

The surface charge density produces the electric field of 
its own, which polarizes the dielectric to form an inhomo- 
geneous polarization P(r) near the interface. It decays to 
the uniform polarization field P, associated with a uni- 
form external field E cx t, far from the interface. Since the 
spatial extent of the interface is neglected, the polariza- 
tion P(r) extends continuously up to the mathematical 
dividing surface. If the dielectric borders vacuum, the 
polarization field changes discontinuously from P(r$) at 
the dielectric side of the surface to zero at its vacuum 
side. 

The sum of the electric field from the surface charge 
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density and the external field is the Maxwell field E(r). 
The Maxwell field, as well as the polarization P(r), can- 
not be directly measured, but can be retrieved from elec- 
tric fields inside cavities carved in the dielectric, as orig- 
inally suggested by Thompson and Maxwell.—^ For the 
simplest geometry of an empty spherical cavity, the so- 
lution of the Laplace equation with Maxwell's bound- 
ary conditions gives the field at the cavity center (cavity 
field^ 

= 2 .^ ext; (1) 

where e is the dielectric constant of the medium sur- 
rounding the cavity. The qualitative prediction of this 
result is that external field E cxt is diminished by a factor 
of 3/ (2e) inside cavities created in highly polar dielectrics 
such as water. This scenario then describes a well-defined 
physical setup testable by laboratory or numerical exper- 
iment. 

One wonders to what extent the mathematical formal- 
ism of Maxwell's electrostatics applies to interfaces of po- 
lar liquids. The interface is obviously not a mathematical 
surface, but has a finite width (Laplace vs Poisson defi- 
nition of the interface^) . In addition, liquid dipoles have 
the freedom to rotate and adjust to the lack of molecu- 
lar interactions from the cavity side of the interface. For 
polar liquids, this restructuring results in preferential in- 
plane orientation of the liquid dipoles at free planar and 
closed interfaces Unless the external field orients the 
dipoles off-plane, such orientational structure eliminates 
the surface charge since op = P n ~ (Fig. [TJd). The 
standard boundary conditions of Maxwell's electrostatics 
do not apply, thus affecting the observable cavity field. 

Water presents a particularly important test case for 
understanding the interfacial electrostatics of polar liq- 
uids. Strong hydrogen bonds between surface waters can 




FIG. 1. Cartoon of the charge distribution op at the surface 
of a spherical cavity carved from a liquid dielectric. Panel 
(a) represents the picture of standard Maxwell's electrostat- 
ics of dielectrics, when the normal projection of the dipolar 
polarization P along the external field E oxt creates the pos- 
itive and negative lobes of the surface charge distribution. 
Those result from the mathematical surface cutting through 
polarized dipoles of the liquid shown by arrows. The overall 
electric field inside the cavity is then reduced from the ex- 
ternal field by the opposing field of the surface charges [Eq. 
0]. Panel (b) shows the scenario suggested by the orienta- 
tional dipolar order at the surface of a "Rossky cavity" in a 
polar liquid.—^ In-plane orientations of the surface dipoles, 
unaltered by a weak external field, do not produce surface 
charge. The result is zero surface charge density op and the 
field at the cavity center following the Lorentz equation [Eq. 
@], Both sets of lines have been produced with e = 72.2 for 
SPC/E water at 300 K. 



potentially prevent their reorientation not only by a weak 
external field of the dielectric experiment, but also by in- 
ternal fields of solute charges. Since it is the response 
of the water dipoles that determines the free energy of 
hydration^ this problem goes beyond the question of 
measuring the electric field inside a dielectric cavity in 
Maxwell's "gedanken experiment" . 

We have previously approached the problem of liquid 
interfacial electrostatics by literally following Maxwell's 
recipe and measuring, by numerical simulations, the elec- 
tric field inside an empty hard-sphere cavity carved in the 
model fluid of dipolar hard spheresi 14 ' 15 Theses studies 
have shown that the cavity field indeed does not con- 
tinuously decrease with increasing e, as Eq. ([1]) would 
suggest, but instead levels off, as a function of e, at a 
value close to the result well established in the theory of 
dielectrics. This limit corresponds to the Lorentz virtual 
cavity^ The latter is defined as a part of the dielectric 
confined by a mathematical closed surface inside it. The 
Lorentz cavity field is then the field produced by the di- 
electric outside this surface. Since no physical interface is 
present, there is no physical polarization at the interface 
and up = by definition^ 

The electric field at the center of a virtual cavity is 
obtained by integrating over the field contributions from 
homogeneous polarization of the medium, instead of in- 



homogeneous polarization in the case of Maxwell's elec- 
trostatics (see below). The field accumulated by this uni- 
form polarization results in the Lorentz field at the cavity 
center 

El = ' E cxt - (2) 

The main qualitative difference between this result and 
the standard cavity field given by Eq. (J) is that the 
Lorentz field does not decay to zero at e — > oo and instead 
levels off at a non-zero value of Eh/E cxt — > 1/3. 

The model dipolar system studied previously by u a 14 i 15 
has in principle shown that interfacial orientational struc- 
ture of a polar liquid might be closer to the limiting sce- 
nario sketched in Fig. |T|b than to the standard, solid-like 
picture sketched in Fig. QJi. The question of practical- 
ity of this observation remains however open. A vacuum 
hard-sphere cavity is a purely theoretical construct. Nev- 
ertheless, the scenario arising from this model can poten- 
tially describe non-polar particles solvated in real polar 
liquids. 

In order to approach this more realistic situation, we 
study here the scenario that we have dubbed the "Rossky 
cavity". The actual system is a Lennard- Jones (LJ) so- 
lute inserted in SPC/E water. As was originally shown 
by Rossky and co-workers^ and supported by many sub- 
sequent studies^ - — water dipoles orient in-plane at in- 
terfaces with non-polar solutes. Although a non-polar 
solute studied here is certainly not a cavity, the orienta- 
tional water structure is highly resilient to external per- 
turbations and remains nearly intact for a broad range 
of solute-solvent attractions^ The term "cavity" is then 
used to stress that it is a large energy of the water hy- 
drogen bonds (AH ~ 4 - 5 k B T per bond for SPC/E 
water—) that supports the interfacial orientational struc- 
ture. Therefore, in contrast to the virtual Lorentz cavity, 
this physical cavity, like the traditional Maxwell cavity, is 
meant to represent realistic measurements of local fields 
inside dielectrics, in water in the present study. The 
Rossky cavity is meant to represent physical situations 
when the interfacial structure is dominated by in-plane 
orientations of water dipoles. We show that the bound- 
ary conditions imposed by molecular interfacial order of 
this cavity produce the electrostatic response approach- 
ing the conditions expected for the virtual Lorentz cavity 
and thus dramatically deviating from the Maxwell cavity 
field. 



II. RESULTS 

The model considered here consists of a single non- 
polar solute interacting with a large number of waters 
mimicking a typical solvation experiment. The solute- 
water interaction is modeled by a Kihara potential com- 
bining a hard-sphere core of the radius i?us with a surface 
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FIG. 2. A cartoon of the Kihara solute in SPC/E water and 
the definition of the projection angles describing the orienta- 
tions of waters in the first hydration layer in respect to the 
surface normal f = n. 9b is the angle between O-H bond and 
f , 9m is the angle between the water dipole and f , 9± defines 
the orientation of the plane of H2O (note that rh • f ± — 0). 



LJ layer of width cqs (Fig. [2]) 
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Here, "0" and "s" are used to label the solute and solvent 
(water), respectively; eo s is the energy of solute-solvent 
LJ attraction. Further, in order to map hard-sphere cavi- 
ties used in our previous studie s 14 ' 15 on the "soft" Kihara 
solute, we adopt the radius of the closest approach 



Rqs — Ras + °os 



(4) 



where tros = 3.0 A has been used in all numerical calcu- 
lations. 

The Kihara potential is the only solute-solvent interac- 
tion introduced in our present model. We are therefore 
missing the effect of electronic polarizability of the so- 
lute always present even for non-polar particles. A crude 
way to estimate this effect is to notice that only the ra- 
tio of the dielectric constants of two interfacing media 
matters for the interfacial dielectric response in dielectric 
theories^ While this prescription requires testing for mi- 
croscopic liquid interfaces, the model situation of a fully 
non-polarizable solute can be mapped on real situations 
of solvated non-polar particles by rescaling the solvent 
dielectric constant with the dielectric constant eo of the 
solute, e — > e/eo, if such property can be reasonably de- 
fined. We will therefore proceed with the present model 
using the term "cavity field" to describe the electric field 
at the center of the solute. 

Configurations of the water-solute mixture were pro- 
duced by Molecular Dynamics (MD) simulations of a sin- 
gle solute inserted in a box of SPC/E waters at standard 
conditions (Supplementary Material (SM) 1 ^). The ratio 
of the cavity and external fields is obtained in the linear 
response approximation as the correlator of the electric 



FIG. 3. The cavity field inside Kihara solutes (filled points, 
connected by the dashed lines) compared to hard-sphere cavi- 
ties inside dipolar fluids (open points, connected by the dash- 
dotted lines).— Results of MD simulations in the main panel 
refer to two values of the solute-solvent LJ attraction, eos = 
0.65 kj/mol (circles) and 20.0 kj/mol (squares). The results 
for the Kihara solutes are plotted against the radius of clos- 
est approach Rq s /(t s defined by Eq. (a s — 2.87 A). Open 
points in the main panel and the inset refer to sets of data 
obtained for varied reduced dipole moment (m*) 2 = f3m 2 /a^ 
of the dipolar hard-sphere fluids: 0.5 (diamonds), L0 (trian- 
gles), and 3.0 (squares). The results for hard-sphere cavities 
in dipolar fluids are plotted against the radius of the closest 
hard-sphere solute-solvent approach Ro s /a s ; a s is the diam- 
eter of the solvent hard spheres. The dashed horizontal lines 
in the main panel and in the inset refer to the Lorentz result 
(L) of Eq. ((21). The dash-dotted horizontal line in the inset 
refers to the Maxwell result (M) of Eq. Q. 



field E s produced by water at the solute center with the 
water dipole moment M s 



E c /E mt = 1 + (0/3) (5E, ■ SM S ) - E 



(5) 



Here, 5~E S and 5M S denote deviations from the corre- 
sponding average values and /3 = l/iksT) is the inverse 
temperature. In addition, the correction term 2S corr in 
the rhs of Eq. ([5]) accounts for the cutoff of the electro- 
static interactions specific to a given simulation proto- 
col. The result for Ewald sums is given in the SM^ 
while the reaction field protocol was covered in previous 
studies, 1 ^ 

Equation ([5]) was used to calculate the field inside the 
solutes of varying size i?ns of the Kihara hard-sphere 
core, and the results are plotted in Fig. [3l In order to 
asses the effect of a uniform solute-solvent attraction on 
the cavity field we have simulated the configurations at 
two values of the LJ attraction, eo s = 0.65 kj/mol equal 
to the LJ energy between oxygens of SPC/E water and 
eo s = 20 kJ/mol close to the energy required to break 
two hydrogen bonds in bulk SPC/E water.— 

The smaller value of the LJ attraction perturbs little 
the water structure at small values of Rhs > but results in 
a weakly dewetted interface^ at the end of the scale of 
i?HS values. This is meant to imply that the first peak of 
the solute-solvent radial distribution function falls below 
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FIG. 4. The solute-solvent distribution functions gjfj{r) [Eq. 
([B]|] vs the distance from the surface of the Kihara solute 
(R 0s = 13 A). Shown are the solute-oxygen (solid lines) and 
solute-hydrogen (dotted lines) radial distribution functions 
(£ = 0) and orientational distributions for I = 1 (dash-dotted 
lines) and 1 = 2 (dashed lines); (a) refers to eo s = 0.65 kJ/mol 
and (b) refers to eos = 20.0 kJ/mol. 



the corresponding peak of the solvent-solvent distribution 
function. On the contrary, the larger value of the solute- 
solvent LJ attraction produced a substantial increase of 
the surface density of water as judged from the solute- 
solvent pair distribution function (Fig. [4]). Furthermore, 
the radial distribution function is nearly zero between 
the first and second hydration shells indicating layering 
of water at the interfaced 

In order to characterize the orientational structure of 
the interfacial water dipoles, we define a series of distri- 
bution functions recognizing the radial symmetry of the 
problem and producing increasing symmetry orders of 
the water dipoles rhj in projection on the outward nor- 
mal to the spherical solute surface f j at the position rj 
of a water molecule (oxygen coordinates are used for the 
center of mass). The radial distribution functions are 
then defined in terms of border Legendre polynomials 



Pi(cos9 mj ), cos 



irio • f , as follows 
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Here, V and iV are the system volume and the number 
of particles, respectively. The zcroth-order distribution 
function is then the standard solute-oxygen radial dis- 
tribution 5qs ( r ) = 9os(r)- In addition to these radial 
distributions shown in Fig. [H we also calculate the or- 
der parameters of the dipoles in the first hydration layer 
by integrating the radial functions over the volume V 1 
defined by the condition < r < Rq s + a s /2 
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Here, p s is the water number density and the integrated 
radial function is normalized to the number of waters in 
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FIG. 5. The first (upper panel) and second (lower panel) 
orientational order parameters p{ 2 of the first-shell SPC/E 
waters vs Rq s = -Rhs + o"0s- The solid circles (eos = 0.65 
kJ/mol) and squares (eos = 20 kJ/mol) refer to Kihara solutes 
in water at T — 300 K. The corresponding open points refer 
to T = 273 K. The crosses in the lower panel show p\ for 
HS cavities in the fluid of dipolar hard spheres (DHS) 22 with 
the reduced dipole moment {m*) 2 = /3m 2 /crj = 3.0; m is the 
dipole moment and a 3 is the HS diameter of the solvent. The 
filled diamond labeled "plane" marks pi for a planar water 
interface from Ref. 0. 



the first hydration layer 
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The results for first-order, p{, and second-order, p^, order 
parameters vs the solute size are plotted in Fig. [5] 

The structure of surface waters is a "squashed" hexag- 
onal ice lattice in which two oxygen sublayers of hexag- 
onal sheets arc brought into one plane of the closest 
solute-solvent approach corresponding to the first peak 
of the solute-oxygen radial distribution function (Fig. 
|4]). The hydrogens of first-shell waters are randomly dis- 
tributed for weak solute-solvent attraction of a hydropho- 
bic surface^ (Fig. [4^,) , but get ordered with increasing the 
attractive pull of the solute (Fig. [4})). The breaking of 
the sublayers of the hexagonal ice sheets results in buck- 
ling of the O-H-0 bond32 from the straight line found 
in hexagonal ice to preferential angles of ~ 10° (more 
populated) and ~ 70° (less populated) (see the SMi2). 
This distribution of hydrogen-bond angles is also quite 
insensitive to the strength of the LJ attraction. Further, 
consistent with a broad distribution of the first-shell hy- 
drogens shown by the dotted lines in Fig. 01 the water 
dipoles are broadly distributed for weak solute-solvent 
attraction, but become more ordered when LJ attrac- 
tion becomes stronger (Fig. [6)3). In all cases, however, 
the in-plane orientation of the water dipoles is preferred, 
and is highly resilient to the changes in the solute-solvent 
attraction. 
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FIG. 6. Distributions of three angles used to define the water 
orientations relative to the radial direction (Fig. [2]): (a) 8t for 
the O-H bond (cos 6), — fj • r), (b) 9 m for the water dipole 
(cos0 m = rh ■ f), and (c) 9± for the normal to H2O plane 
(cos#^ = fj_ ■ f); R Qs = 13 A. 



The distributions of projections of water O-H bonds 
and the normal to the H2O plane (Fig. [2]) are shown in 
Fig. [6K,c. Both are consistent with the picture of the 
in-plane orientation of the water dipoles, with the first- 
shell waters populating variably the states of H2O plane 
perpendicular to the surface normal and slightly tilted, 
at ~ 10° in respect to the normal. These two states are 
populated in the ratio of about 1:2 for the weak attrac- 
tion, but the tilted state becomes a dominant one at the 
stronger attraction. This observation is consistent with 
the general tendency of interfacial waters to increasingly 
occupy one particular orientational state with increasing 
strength of attraction of either protons or oxygens to the 
solute. In contrast, the hydrogen disorder generally in- 
creases when a combination of two types of attractions 
are present at the interface^ a situation typical of hy- 
drated proteins. 



III. DISCUSSION 

The results for the electric field inside the Rossky cav- 
ity in water are shown by filled points in Fig. [3] They 
clearly tend to the Lorentz field limit. We also plot in 
Fig. El (open points) the results for the electric field inside 
hard-sphere cavities in dipolar fluids . 14 ' 15 Although the 
two situations are physically distinct, a non-polar solute 
in the former case vs an actual cavity in the latter, the 
phenomenology of the polar interfacial response is generic 



for both cases. The common output of these two sets 
of simulations suggests that the picture of non-polarized 
interface represents the response of free and non-polar 
interfaces of polar liquids in general. 

It is instructive to see how the Maxwell cavity field 
and the Lorentz field appear based on the assumptions 
regarding the surface charge density (Fig. [TJi,b) . We first 
assume a constant value for the surface charge density 
op = <Jq. The Maxwell electric field in the medium is 
then the sum of the external field and the radial field 
propagating from a uniformly charged interface 



E(r) 



1 47Tgggo 

--Ejcxt H , — r. 



(9) 



This electric field will polarize the dielectric- yielding the 
polarization field P(r) = (e — l)E(r)/(47r). This polar- 
ization field is inhomogeneous close to the interface and 
decays to the uniform polarization P in the bulk. The 
polarization of the dielectric will in turn produce its own 
electric field, which, combined with the external field, 
gives the cavity field 



E r — E r 



/ T(r) • P(r)dr. 

hi 



(10) 



Here, T(r) = VVr" 1 is the dipolar interaction tensor 
and the integral is over the dielectric occupying the vol- 
ume Q outside the dielectric cavity of radius Rq. The 
radius Rq does not need to be specified since the calcu- 
lation results do not depend on its value. 

The radial polarization arising from the second sum- 
mand in Eq. ([9]) gives zero contribution to the cavity 
field and the result of Eqs. ([9]) and (fT0|) is the Lorentz 
field given by Eq. @. A constant, angular- independent 
surface charge density makes therefore no contribution 
to the cavity field. It however contributes to a non- 
zero, spatially constant electrostatic potential inside the 
cavity^ 5 - 

Given axial symmetry of the problem, surface charge 
density can be expanded in Legendre polynomials, 
ap(9) = J2e GtPi (cos 9), where 9 is the polar angle with 
the direction of the external field. Since the zero-order 
term ctq does not contribute to the cavity field, one can 
take the first-order term and, neglecting the quadrupo- 
lar and higher moments, have the dipolar approximation 
<jp(9) = o\ cos 9. The dipole moment created by the so- 
lute interface is then M = <7if2o • This is the situation 
sketched in Fig. QJi. The direct solution of the Laplace 
equation with Maxwell's boundary conditions results in 
the interfacial dipole 



Mn = -P 



2e + 1' 



(11) 



where P = (e — l)/(47re)E cxt . The negative sign here 
indicates that the interface dipole Mo orients oppositely 
to the external field, as is also clear from Fig. QJi. 

The dipolar surface charge density following from Eq. 
(fTTj) is shown by the solid line in Fig. This surface 
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FIG. 7. Surface charge density ap(9) — P n (8) vs the polar 
angle 9 measured from the direction of the external field (a). 
The solid line (M) shows the result of Maxwell's electrostatics 
[Eq. (H6|), e = 72.2 for SPC/E water at 300 K]. The dots refer 
to up calculated as a linear response to an external field from 
40 ns of MD trajectories (see the SM±2). The MD a P is 
collected from a water layer Ros < r < Ro s + 0.1a s for the 
solute with the size Ros = 13 A and eos = 0.65 kj/mol. The 
dashed line shows the result of using E c from simulations to 
calculate xi from Eq. (|13[) . Panel (b) shows the profile of 
water's dielectric constant e(r) = 1 + (4TrP/3Q(r)){S~M. a (r) ■ 
8M S ) (see text). 



charge density creates a dipolar field, which adds to the 
external field to produce the Maxwell field 
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When the polarization field P(r) is calculated from Eq. 
(fl~2j) and substituted into Eq. (fT0)l . the cavity field be- 
comes 
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is the dipolar response function of the solute interface 
to the polarizing external field. The Maxwell result for 
M in Eq. (fTTj) gives the standard cavity field of Eq. ((T|) , 
while xi = M = leads to the Lorentz field of Eq. J2]). 

Surface charge density ap(9) provides a convenient 
mathematical idealization of the polarization of the in- 
terface. Its direct calculation from real-space, finite-size 
simulations presents, however, a significant challenge. 
The problem is illustrated in Fig. [7K which shows op{9) 
calculated at the solute surface from MD trajectories us- 
ing the linear response approximation (see the SMi^). 
The result is a clearly angular-independent function, pro- 
ducing o"x — and E c ~ E\^. This outcome is, however, 
not entirely consistent with a slightly positive o\ follow- 
ing from the substitution of the simulated E c into Eq. 
(|T3)) (dashed line in Fig. [3i). 

The reason of possible uncertainties in the attempts 
to calculate ap{9) from sampling surface dipole direc- 
tions from MD trajectories is illustrated in Fig. (7}d. The 



definition of a mathematical dividing surface is uncer- 
tain when the dipolar response is strongly varied at the 
interface, as is the case here. We show the dielectric con- 
stant of the polar layer of radius r from the cavity cen- 
ter obtained from the linear response approximation as2£ 
e(r) = l + (4itP/3Cl(r)){SM s (r)-dM s ), where M s (r) and 
51(r) are, respectively, the dipole moment and volume 
of the solvent within the r-surfacc. The layer dielectric 
constant peaks near the interface and then slowly de- 
cays to the bulk dielectric constant e. Any surface drawn 
next to the interface will therefore reflect a different vari- 
ance of the surface dipole, with a different outcome for 
(jp{9). How much the functional form of ap(9) changes 
depending on the surface definition cannot be answered 
here. Attempts to obtain ap(9) at surfaces deeper into 
the bulk did not produce converged functions and there- 
fore are not shown here. 

It seems worth emphasizing some critical differences 
between both (Maxwell and Lorentz) continuum results 
and the molecular picture offered by the present numeri- 
cal simulations. The continuum cavity field does not de- 
pend on the cavity size, but only on the cavity shape (e.g., 
on the length/diameter ratio for a cylindrical cavity^). 
This is of course the consequence of neglecting the ac- 
tual spread of the interfacial region relative to the size 
of the solute, which is behind the definition of the inter- 
face as a mathematical dividing surface. Our simulations 
allow us to estimate the solute size at which this approx- 
imation becomes valid. Figure [3] shows that the cav- 
ity field starts leveling off to a size-independent limit at 
Ro s /a s ~ 2 — 2.5. It gives an estimate of the solute size 
at which the "Lorentz continuum" starts to take hold: 
the solute needs to be 4 — 5 times larger than the solvent 
molecule. The "Maxwell continuum" is, however, never 
reached in our simulations. 



IV. EXPERIMENTAL OBSERVABLES 

The surface polarization effects discussed here have a 
number of consequences observable in a macroscopic lab- 
oratory experiment. The total free energy of polarizing 
the dielectric is measured by the dielectric experiment 
and is affected by the dipole accumulated at the inter- 
face with a non-polar solute. The total dipole of a mix- 
ture sample -M m ; x is reduced relative to the homogeneous 
solvent by the volume excluded by the solute and, in ad- 
dition, is affected by the dipole of the interface 



M tl 



= PO-(2/3)(e-l)M iVo 



(15) 



Here, No is the number of the solutes in the solution 
and fi, as above, is the solvent volume. The solutes are 
assumed to be non-interacting, although the theory can 
be extended to non-ideal solutions requiring the solute- 
solute structure factor as an additional input<21 

The standard arguments of the theory of dielectrics 
then suggest the equation for the dielectric constant of 
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FIG. 8. Relative dielectric constant increment Ae/e, Ae = 
£mix — e for several aqueous solutions: lysozyme 29 (filled 
squares), dioxane— (open circles), glucose^ (open diamonds), 
and maltotriose^ (open triangles). The solid line (M) shows 
the result of Eq. (|18l) . while the dashed line (L) is the di- 
electric increment for a solute with no surface charge density, 
o\ = 0, and therefore Lorentz result for the cavity field [Eq. 
(|13|) ]. Experimental and calculation results used to convert 
the experimentally reported solute concentrations to volume 
fractions can be found in the SM.— The dash-dotted lines 
connect the experimental points for two saccharides. 
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6mix 



l + Tto(e-l)(l + (87r/3)exi), 



(16) 



where % = N n /(N £lo + ^) is the volume fraction 
of the solutes. Equations (fl3|) and (fT5)l . taken together, 
suggest that dielectric constants of low-concentration so- 
lutions and cavity fields inside solutes both give experi- 
mental routes to measure the dipolar response function 
of the interface xi- Combining two equations together, 
one gets 



= l + 377o 



- 1 



(17) 



For the cavity field E c = Eyi given by Maxwell's form 
in Eq. ([T]) one arrives at the resul t 15 ' 27 



1 



3(6-1) 
2e+l ' 



(18) 



This equation is consistent with the Maxwell- Wagner 
theory 2 ^ in the limit of small volume fraction t]q . 

Figure [8] illustrates the application of Eqs. (fT6|) and 
(|18| to solutions of ionic and hydrogen-bonding sub- 
stances in water^2r— The closed squares in the figure 
show the relative dielectric constant increment Ae/e, 
Ae = e m jx — e for hydrated lysozyme^ At the pH of 
the measurements the protein carries the total charge of 
+10 electron units. Open circles in the figure show the 
results for the dioxane-water mixture. 30 Dioxane offers 
strong hydrogen-bond acceptor sites. Finally, diamonds 
and triangles show, respectively, hydrated glucose 3 ^ and 
maltotriose ^ both offering multiple hydrogen-bonding 



sites. The results for lysozyme and for two saccharides 
refer to dielectric increments of the loss peak of the wa- 
ter component of the solution ! 29 ' 31 This approach has al- 
lowed us to eliminate the contribution of the permanent 
dipole moment of the solute, characterized by a much 
longer relaxation time. The increment of the water loss 
peak reflects only the interfacial polarization described 
by the present model. The dipole moment of dioxane is 
zero and the increment of the static dielectric constant 
was used in that case. 

The experimental results are compared to Maxwell 
[Eq. (fpgj) . solid line] and Lorentz [xi = in Eq. (fToj) . 
dashed line] scenarios. All solutes, potentially breaking 
the interfacial network of hydrogen bonds, arc cither con- 
sistent with the Maxwell interfacial polarization or fall 
between the Maxwell and Lorentz predictions. 

The Lorentz interface predicts a much enhanced sensi- 
tivity of the dielectric constant to non-polar impurities, 
which also implies a significant free energy penalty for 
polarizing such mixtures. Experimental testing of the 
Lorentz scenario should therefore meet with obvious sol- 
ubility difficulties. However, scenarios intermediate be- 
tween the Lorentz and Maxwell limits clearly exist and 
the concept of an effective surface charge density up pro- 
vides a useful conceptual framework for developing theo- 
ries of polar response not restricted to Maxwell's bound- 
ary conditions. 

While uniform fields of the dielectric experiment mea- 
sure the dipolar response of the entire solution, non- 
uniform external fields give access to the response func- 
tion xi an d therefore to the dipolar polarization of the 
interface. If the external field varies on a scale large com- 
pared to the dimension of the solute, the interfacial dipole 
couples to the field gradient. The result is a force acting 
on the solute 



F z = (f2 /2)xiV z £: ox t(z) 2 



(19) 



where the z-axis is chosen along the external field. This 
phenomenon, known as dielectrophoresis^ allows a di- 
rect access to the polarization of the interface. The 
Maxwell solution for xi [Eqs. (fTT|) and (fTf]) ] leads to the 
standard dielectrophoresis coefficient K oc xi used in the 
theory of colloidal suspensions^ 2 - In these applications, e 
is typically replaced with the ratio of the dielectric con- 
stants of the solvent and the solute, e — )• e/eo, with the 
result 



eo 



e + 2e 



(20) 



If the solute is more polar than the solvent, one gets 
positive dielectrophoresis (K > 0, xi > 0), an d solute's 
attraction to a stronger electric field. The opposite case 
of K < (xi < 0) implies negative dielectrophoresis and 
thus repulsion of the solute from the field. 
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V. SUMMARY 

This paper extends our previous studie d 14 ' 15 of dipolar 
fluids interfacing hard-sphere cavities to attractive non- 
polar solutes hydrated by SPC/E water. In-plane orien- 
tational structure of the surface dipoles holds for inter- 
faces of both dipolar liquids and water. The oricntational 
distribution of the intcrfacial dipoles alters the boundary 
conditions of the polar response problem monitored here 
in terms of the electric field inside the solute. Like in 
previous studies of cavities in dipolar fluids, the field in- 
side the solute does not follow the prediction of Maxwell's 
electrostatics [Eq. (JTJ)] and instead tends, with increasing 
solute size, to the limit established by the Lorentz field 
of a non-polarized interface [Eq. J2|)]. 

We find that the deviation of the cavity field from the 
Lorentz result is given by the dipolar response function 
of the interface that also enters the dielectric constant 
of an ideal solution. Finding the slope of the dielectric 
increment of the low-concentration mixture vs the solute 
concentration [Eqs. (fT7|) and (|16p ] thus provides a direct 
input into the cavity field [Eqs. (|17[) and (flU)) ]. This state- 
ment also applies to the frequency-dependent response. 
As such, the dielectric increment of the fast water com- 
ponent of the solution was taken in Fig. [8] to avoid the 
effect of a slowly relaxing permanent dipole of the solute. 
Our formalism is therefore applicable to ionic and polar 
solutes when the frequency range exceeds the loss peak 
of the solute dipole. Along these lines, the formalism can 
also be used at much higher frequencies of UV/VIS light 
to find the refractive index corrections for rates of radia- 
tive decay of photoexcited chromophores^ and quantum 
dots^ in solution. According to our present results, this 
input should be sought from measuring the refractive in- 
dex of corresponding mixtures [Eq. (|17j) ]. 

We also find that small differences in the cavity field, 
which arc hard to resolve within the current simulation 
protocol, propagate in very substantial differences in the 
dielectric response of a mixture. While this fact also 
spells out a thermodynamic difficulty of preparing solu- 
tions that might test the Lorentz interface by conven- 
tional dielectric spectroscopy, a large difference in the 
slopes of dashed (Lorentz) and solid (Maxwell) lines in 
Fig. [8] leaves much space to looking for intermediate sce- 
narios. 

Returning to the question posed in the Introduction, 
the liquid interface assumes the structure that makes its 
contribution to the field inside and outside the solute es- 
sentially null. The whole solvent response to an external 
polarizing field is given by the uniform polarization of 
the bulk and does not include an inhomogeneous com- 
ponent due to the surface charge [second summand in 
Eq. (|12p]. In retrospect, this outcome should have been 
expected. The standard practice of liquid state theories 
and corresponding numerical simulations suggests a short 
propagation length of perturbations in liquids. A cavity 
or a non-polar solute should therefore be "invisible" to 
an observable (e.g., Maxwell field) measured a few sol- 



vent diameters from the interface, as indeed our results 
show. This result of course goes against the concept of 
the surface charge inducing a long-ranged Coulomb per- 
turbation in the solvent. 

The original concept of electric polarization envisioned 
by Maxwell anticipates a limited, clastic response of 
medium's electric charge to an external electric field.- 
For molecular dipoles of a dielectric material, this con- 
cept implies limited small-amplitude reorientations of the 
dipoles, aligning them with a weak external field (linear 
response). This physical picture, which seems to match 
the problem of dipolar polarization of a free solid-like 
interface, was historically extended to liquid dielectrics 
which, in contrast, have more extended ability to re- 
spond by both changing the positions of their dipoles 
and producing large- amplitude dipolar reorientations. In 
other words, the translational and rotational mobility 
supported by the liquid phase allows the surface dipoles 
to react to the creation of the interface by rearranging 
their dipolar orientations in a way of diminishing the 
stress of a sharply varying density profile. This orienta- 
tional order responds to a weak external field by the rules 
that do not require a surface charge and the induction of 
an intcrfacial dipole combining the negative and positive 
lobes of the surface charge density (Fig.QJi). The bound- 
ary conditions applicable to these systems differ from the 
ones anticipated for Maxwell dielectrics. 

By the way of yet another historical aberration, the 
standard theory of dielectrics has been mostly applied 
not to free liquid interfaces, for which Maxwell's con- 
struct was put forward, but to highly hydrophilic and 
wetted interfaces covered with surface ions or to solvation 
of molecular ions (Born theory of ion solvation and its ex- 
tensions). The agreement between observations done for 
these surfaces and solutions with the predictions of di- 
electric models is often used to support the conceptual 
framework of dielectric theories developed for free or non- 
wetted (small attraction) surfaces. Our present develop- 
ment cautions against this inconsistency and points out 
that the orientational structure of the interface defines 
boundary conditions of dielectric theories and with that 
the microscopic and macroscopic fields observed near di- 
electric interfaces. 

The picture of in-plane dipolar polarization of the liq- 
uid interface is rather robust and insensitive to the inter- 
face curvature. It holds even for planar surfaces (filled 
diamond in Fig. [S]^^ 5 - suggesting that the conclusions 
reached here for an admittedly narrow range of solute 
radii may extend to larger solutes of potentially meso-to- 
macroscopic size. 

What matters for the boundary conditions entering the 
electrostatic response functions is the orientational distri- 
bution of surface waters. This can be altered by surface 
ions and polar group o 12 ' 36 " — which potentially can cre- 
ate a non-zero surface charge density up ^ matching 
the standard conditions of Maxwell's electrostatics. The 
standard prescriptions will apply to those cases. It fol- 
lows directly from Coulomb's law that the change of the 
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normal component of the electric field at the interface of 
a dielectric with vacuum is related to the surface charge 
density, —AE n = 4nap. The near-zero up then implies 
the continuation of the normal component of the electric 
field, and, therefore, the continuation of the vector E of 
the Maxwell field across a dielectric interface. 

Liquids with large cohesive energy, network liquids, 
such as water, in particular, seem to be particularly 
relevant to this discussion. The strength of hydrogen 
bonds (^4 — 5 k^T per boncU^) is so significant that in- 
plane dipolar pattern may withstand local electric fields 
of solute partial charges. The ratio of the characteristic 
strengths of the solute-solvent to solvent-solvent interac- 
tions will therefore determine the orientational distribu- 
tion of the surface dipoles and, ultimately, the type of 
boundary conditions used in calculations of the electro- 
static response. 
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